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Abstract 

The current theoretical status of the Lamb shift in He"*" is discussed. Recent calculations of two- 
loop binding corrections to the Lamb shift significantly shift the theoretical value of the "classic" 
Lamb shift in He^, i.e. of the 2Si/2~2Pi/2-iiiterval. In this brief research note, we present a new 
(theoretical) value for this interval which reads 14041.474(42) Mhz. The theoretical uncertainty 
is reduced as well as the discrepancy between theory and experiment. Planned measurements 
should be of help in further elucidating the situation. 

1 Introduction 

One of the intriguing questions which remain to be answered in bound-state quantum electrodynamics 
is related to the discrepancy between the theoretical and experimental values for the Lamb shift in 
ionized helium, or He"*", which is a hydrogenlike atomic system with a nuclear charge number Z = 2. 
The latest measurement of the "classical" Lamb shift in He"*", using the anisotropy method, has led 
to the result 

i^exp(2Si/2 - 2P1/2) = 14041.13(17) MHz . (1) 
The theoretical value given in |2] is different: 

i^th,oid(2Si/2 - 2P1/2) = 14041.57(8) MHz , (2) 

where the uncertainty has been assigned according to an estimate of the magnitude of higher-order 
contributions which had been unknown. If these higher unknown terms had turned out to increase 
rather than decrease the value of the Lamb shift, then a 3.0 a discrepancy between theory and experi- 
ment could have arisen. In the current paper, we review a shift in the theoretical value due to recently 
calculated higher-order two-loop self-energy corrections, and a further shift due to radiative-recoil 
corrections where a certain discrepancy between conflicting results has recently been eliminated. 

Bound-state quantum electrodynamic (QED) energy shifts may be expressed as combined perturbation 
series in the parameter a, which is the fine-structure constant, and in Za, which is a measure of the 
electron-nucleus interaction. Although Lamb-shift measurements in atomic hydrogen (Z = 1) have 
a long history, it is necessary to perform at least one measurement for a different value of Z in 
order to test the Z-dependence of the corrections, which are generated by bound-state effects and are 
fundamentally different from the QED of free electrons. In He+, we have Z = 2, and the series in Za 
converge much slower than for Z = I (atomic hydrogen) . 

A more accurate measurement of the Lamb shift in He+ is currently being pursued by Hinds and 
Boshier It is perhaps worthwhile to note that He+, as it constitutes an ionic system, should also 



be an attractive candidate for high-precision laser-spectroscopic measurements in traps, where the 
small Doppler shifts associated with the slow movement of the trapped ion could provide a basis for 
a further significant reduction in the experimental uncertainty. 



2 Evaluation of Higher-Order Corrections 

A detailed summary of radiative corrections which contribute to the Lamb shift of n = 2 states in 
He"'" has been given in Table I ofpp. Recently, the understanding of two-loop binding corrections to 
the Lamb shift of S states has been significantly enhanced |21 ■ Terms of order 

" '^-^^^^-^ln\Za)--]B,,{nS) and (^^)' ^^^^1!^ InKZa)-] i?ei(nS) (3) 

have been evaluated. [Here, mo is the electron mass, c is the speed of light, and the energy shift can 
be converted to frequencies via division by the natural unit of action, which is Planck's constant h. In 
the following, we use a system of units in which h = = c — 1. We will also use the notations m and 
nic for the electron mass interchangeably Jji, and by A/, we denote the mass of the a particle.] These 
coefficients involve contributions originating from two-loop self-energy, combined self-energy vacuum- 
polarization, and two-loop vacuum-polarization diagrams. The final results Q read as follows, 

S62(1S) - 1^-1^ ln2- -0.461 891, (4a) 
135 9 
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Be2{nS) = B62(1S) + — (- + — -- -ln(n) + *(n) +7) , (4b) 



where 7 — 0.57721... is Euler's constant, and ip{z) — r'(2:)/r(z). For n ^ 00, the coefficient 
approaches a limiting value of 

lim Be2{nS) = -i^ [71 + 60 (7 - In 2)] = 1.897 603 . (5) 

rwoo 135 

The numerical value of the i362-coefficient depends rather significantly on n. For n — 2, we have 

419 32 

Be2 (2S) = — - — In 2 = 0.639 180 . (6) 
loo y 

A graph of BQ2{nS) in the range 71 = 1, . . . , 8 is given in Fig. ^ 
Surprisingly large results have been obtained .2, for the coefficient -Bgii 

ln2 + ^ ln22 + ^iV(lS) = 50.344005, (7a) 
135 9 3 

BeiinS) = i?6i(lS) + ^ [N{nS) - iV(lS)] (7b) 
304 32 , \ /3 1 1 , , , ^, , 

135-yi^'j U + i^"^"'"^"^ + *^"^+^ 

The quantity N{nS) results from a correction to the Bethe logarithm induced by a local potential, 
and has been given for excited states in In particular, we have A^(2S) = 12.032 209(1) and 
i36i(2S) = 42.447 669(1). This is a 15% deviation from the corresponding result for IS which 
is indicated here in equation H7a|) . The i?6i-correction, for 2S, therefore does not quite enhance the 
2Si/2-2Pi/2 interval as much as might be expected from the corresponding result for the IS state 
[equation H7a|l ] . 
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Figure 1: Dependence of the _Bg2-cocfBcient on the principal quantum number 
n for S states. On the abscissa, we have the principal quantum number n, 
and the ordinate axis represents the numerical values of 3^2(1^8) as given in 
equation Q . The smooth curve results from a three-parameter fit with a model 
of the form a + b/n + c/n^ (the quantities a, b, and c are fit parameters). A 
model of this form has been shown to lead to a satisfactory representation of 
the n-dependence of quantum electrodynamic radiative corrections in many 
cases (see ^ and references therein). The coefficient Bq2 changes sign between 
n = 1 and n = 2. 



Additionally, we note that the recently available ^ nonlogarithmic two-loop terms of order 

/aN2 {Zaf mc 

[-) ^3 BeoinS) (8) 

result in a negative energy shift for S states and that the _Bgo-coefficients are surprisingly large, 

BeoilS) = -61.6(3) ±15%, (9a) 
B6o(2S) = -53.2(3) ±15%. (9b) 

Here, the uncertainty is mainly due to an unknown high-energy contribution to this coefficient. 

Recent progress has also been reported with regard to recoil corrections of first order in the mass 
ratio mc/M , where M is the mass of the atomic nucleus. For hydrogenic systems, a nonperturbative 
evaluation to all orders in Za has been performed in and the results are in agreement with the 
known terms of the Za-expansion of this effect. These are as follows: in the order {Za)^ {rric/M) mc c^, 
we have the Salpeter correction. In relative order {Za), we have a correction obtained by Pachucki and 
Grotch [7]. Finally, in relative order {Za)'^ ln^[(ZQ;)~^], we have a further correction recently obtained 
by Pachucki and Karshenboim, and Melnikov and Yelkhovsky The sum of these corrections, 

evaluated for Z = 2, leads to a satisfactory agreement with the numerical data published in ^ for the 
entire range 1 < Z < 100. The nonperturbative (in Za) remainder can conservatively be estimated 
as 3 kHz, which is three times the mag nitude of the term of order (Za)^ \n^[{Za)-^] (see |H10). This 
term is not the main source of the theoretical uncertainty. 

A discrepancy in the analysis of radiative-recoil corrections of order a (Za)^ (in^jM) rric has re- 
cently been resolved. According to the fully analytic calculation ^UJ, the coefficient multiplying this 
term has the value — 1.324 028/n'^ for S states (here, n is the principal quantum number). In J^, a 
slightly different coefficient of — 1.374/n^ had been obtained, whereas in |lll I12L [T5| . a coefficient of 
— 1.988(4) has been given. Assuming the correctness of the most recent evaluation jTO], the result 
for radiative-recoil corrections of order a (Za)^ (m^jM) (? should be taken as —0.014 MHz rather 
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than —0.035 MHz. This concerns the entry termed "a (Za)^ m/M" in Table I of 1 which, in units 
of arricC^, is of relative order (Zotf' m^jM . 

Finally, we note that the uncertainty due to the nonperturbative (in Za) one-loop self-energy re- 
mainder G's,^(Za) has recently been eliminated \i : the numerical result reads G'se(2Si/2, 2a) = 
—30.644 66(5). For the 2S state of He"*", the following energy shifts result from the results reported in 
equations Q-®: 

^62(28) +0.037 MHz (10a) 

(entry termed a {Zaf h\{Za)-'^ in Table I of [T]) , 
B6i(2S) ^ +0.289 MHz (10b) 

(entry termed a {Zaf \n{Zay^ in Table I of P) , 
B6o(2S) -> (-0.043 ±0.006) MHz (10c) 

(entry termed a {Zoif in Table I of P) , 
radiative recoil -0.014 MHz (lOd) 

(entry termed a {Zcif' ra/M in Table I of J) , 
GsE(2Si/2,2a) ^ -10.624 MHz (lOe) 

(entry termed {Zaf Gs^{Za) in Table I of [T]) , 

For completeness, we should also note slight inconsistencies in the notation of the corrections in ^j. 
First, in Table I of P], all corrections are given in units of a rric (? . The Salpeter correction has inad- 
vertently been termed "(Za)^ m/M" whereas in units of a rric c^, it should be termed "Z {Zaf m/M." 
Correspondingly, the correction denoted "(Za)^moC^" should actually be termed '^Z {Zaf m/M" , 
and the correction denoted as "{Zaf \n{Za) m/M" is of course a correction of (relative) order 
''Z{Zaf In(Za) m/M." 

We now discuss the main source of the theoretical uncertainty. There exists an essentially unknown 
three-loop binding correction of order 

(-)'^^^^^C.o(nS), (11) 

where we note that the Cso-coefficient, in analogy to the two-loop binding correction B^q « 
—21.5562(31), might be numerically large. With the estimate IC50I ~ 30 for the unknown term, 
the corresponding additional uncertainty for the 2S state is 0.008 MHz. The nuclear-size contribution 
to the uncertainty is 0.010 MHz [151 IT] (this is assuming a nuclear radius of 1.673(1) fm). There is 
an additional 0.015 MHz uncertainty due to unknown higher-order two-loop corrections for 2P states. 
We obtain the following theoretical error budget for the theoretical value of the Lamb shift in He+: 

• unknown high-energy part of Bqq: 6 kHz. 

• unknown higher-order two-loop effects for P states: 15 kHz. 

• imknown three-loop binding correction C50: 8 kHz. 

• nonperturbative remainder of nuclear recoil: 3 kHz. 

• uncertainty in the nuclear charge radius: 10 kHz. 

The a particle charge radius has been determined to high accuracy in 15 , using spectroscopic ex- 
perimental data from the system ^"-He^"*". However, the accuracy of the resulting charge radius 
of 1.673(1) fm has been questioned (see e.g. ^JEI)- The result obtained in combining scattering 
experiments is 1.674(12) fm ^H] and has a considerably larger error. The above entry in the error 
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budget should therefore rather be interpreted as a lower limit of this problematic contribution. Be- 
cause theoretical errors cannot be expected to follow a normal distribution like experimental errors, 
there exists no universally adopted procedure for the determination of a total theoretical error in a 
situation where several unknown terms contribute to an error budget. In the current work, we choose 
to (conservatively) add the above errors in order to obtain a total theoretical uncertainty of 42 kHz. To 
complete the discussion, we mention that (i) under the assumption of normally distributed theoretical 
errors, the total theoretical uncertainty would be reduced to 20 kHz, and (ii) under the assumption 
of a 0.012 fm uncertainty in the nuclear radius JH], the theoretical error due to the a particle charge 
radius alone would be 126 kHz. 



3 Conclusions 

We have re-evaluated theoretical predictions for the 2S-2Pi/2 Lamb shift in He+. Taking into account 
recently evaluated corrections to the entries in Table I of U (see equation (|10|l ). we obtain the new 
theoretical value of 

1/^(281/2 - 2P1/2) = 14041.474(42)(126)MHz, (12a) 

where the first error is due to the error budget as discussed above, and the second error results due to 
the nuclear size alone, if we assign a larger nuclear charge radius uncertainty of 0.012 fm 18 instead 
of 0.001 fm jni- For the 2P3/2"2Si/2-iiiterval, the new theoretical prediction is 

i/th(2P3/2 - 2S1/2) = 161542.029(42)(126) MHz, (12b) 

with the same error budgets as for 2Si/2 — 2Pi/2- 

There remains a small discrepancy between the "new" theoretical value of 14041.474(42) MHz (equa- 
tion H12a(l ^ and the experimental result of 14041.13(17) MHz given in equation However, this 
discrepancy is less severe than the difference of the experimental value |^ and the "old" theoretical 
value of 14041.57(8) MHz (see equation Q). 

The sum of the experimental {— 0.17 MHz) and the theoretical uncertainty {— 0.04 MHz) is cr = 
0.21 MHz. This should be compared to the difference 5 = 0.34 MHz between the theoretical and 
experimental "expectation values" l)12a|l and The difference S corresponds to 1.6 cr. If we assign 
the larger error of 0.012 fm to the nuclear radius, then the discrepancy shrinks to 1.2 cr. 

Even though a difference of 1.2 cr to 1.6 ct between theory and experiment remains somewhat unsettling, 
it is not large enough to be regarded as statistically significant. Even at this level, it still represents 
a better test of the unexpectedly large ^50 two-loop binding correction to the Lamb shift (-1.339 
MHz) than that provided by the corresponding measurement in hydrogen. In addition, the anisotropy 
method provides a completely independent method of measuring the Lamb shift that is not limited 
by the large level width of the 2P state (s), and it is the only method that is capable of comparing two 
different atomic or ionic species (H and He+, c.f. Ref. jEl) within the same apparatus. 

It may also be beneficial, in the near future, to re-analyze the results of the most recent measurement ^ 
in terms of QED corrections to the 2S~2P transition matrix elements which are currently under 
study j20| . In the current paper, we restrict the discussion to the theoretical predictions for the 
absolute frequency of the 2S-2P transitions. Further clarification of this and related questions will 
certainly benefit from a reduction of the experimental uncertainty of the Lamb-shift measurement to 
±100 kHz or better, as is currently being pursued by Hinds and Boshier Finally, we stress that 
all theoretical predictions for transitions in He^ rely on the single available accurate determination of 
the charge radius of the ct-particle as described in • There is still a lack of accurate independent 
verifications of this determination of the charge radius. In the end, the answer to a number of questions 
related to Lamb-shift measurements in atomic hydrogen and He+ might be as easy as a slight deviation 
of the proton and a particle charge radii from their currently accepted values. However, other, more 
subtle explanations cannot be ruled out at present, either. 
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